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1. Introduction 

^ The outer crust of a neutron star, below densities ofp^4xl0^^g cm~^ consists of 

matter in a state not too far removed from that found in white dwarfs: a lattice of nuclei 
permeated by a relativistic, degenerate electron gas which gives the dominant contribution 
to the pressure of the matter Q. As pressure increases with depth, equilibrium with respect 
to weak interactions drives the nuclei to become more neutron rich. There comes a point 
when the intra-nuclear forces can no longer bind all the neutrons, and neutron drip occurs. 
Above p ^ 4 X 10^^ g cm~^, a new regime is entered in which the nuclear lattice is bathed in 
a fluid of ('dripped') neutrons. These neutrons are delocalized much like conduction band 
electrons in metals. From this density inwards, the equation of state (EOS) is dominated 
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Figure 1. Cartoon representation of the inner neutron star crust. 

by pressure arising from nucleon-nucleon interactions. 

A cartoon representation of the crustal layers below the outer crust is shown in Fig. 1. 
Models predict two distinct layers. (1) The inner crust between densities p ^ 4 x 10^^ 
g cm~^ and p ^ 10^^ g cm~^ is an elastic solid consisting of a lattice of heavy, neu- 
tron rich nuclei surrounded by fluid neutrons, with the nuclei increasing in size and mass 
with density while the inter-nuclear spacing and nuclear proton fraction decrease ||2l. (The 
presence of a background electron gas will be taken as given from now on). The dripped 
neutrons are expected to become superfluid shortly after the neutron star is formed as it 
rapidly cools below temperatures of 10^ — 10^ K ||3l. (2) The mantle [4] between p ^ 10^"^ 
g cm~^ and the crust-core transition density consists of frustrated matter: the competition 
between the nuclear surface energy and the nuclear and lattice Coulomb energies over sim- 
ilar length scales drives the formation of exotic nuclear geometries termed nuclear 'pasta' 
[5, 6 1 which proceed through a canonical sequence of phases: cylindrical (spaghetti) 
slab (lasagna) cylindrical bubble spherical bubble. The latter three ('bubble') phases 
are distinguished by the delocalization of the charged nuclear component of the matter 
(containing the protons) in one or more dimensions, and corresponding localization of the 
charge-neutral nuclear component (fluid neutrons). Similar microscopic structures are 
observed in terrestrial soft condensed matter systems such as surfactants ||7l[8l; by analogy, 
we can expect rich mechanical properties, intermediate between liquid and elastic solid, to 
emerge in the mantle |9|. Some crust models predict the absence of the mantle 1 10|; its 
presence depends sensitively on the nuclear microphysics of the crust. Fig. 1 gives a range 
of widths for the inner crust and mantle taken from the model presented in this chapter and 
encompassing a range of neutron star masses from 1-2 Mq and a range of equations of state 
as discussed later. 

To describe the states of matter in a neutron star, one needs a model for the nucleon- 
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nucleon interactions as they are manifested in a many-nucleon context. A useful concept 
that bridges the gap between ab initio nucleon-nucleon calculations, nuclear experimental 
observables, and neutron star matter is that of uniform nuclear matter (NM). This is an 
idealized system, homogeneous and infinite in extent, of neutrons and protons interacting 
solely via the strong force. The energy per particle of such a system at a density p and 
proton fraction x, x), is referred to as the nuclear matter equation of state (NM EOS). 
In the regions of the neutron star core where protons and neutrons exist, the NM EOS 
can be combined with the electron energy, and under conditions of charge neutrality and 
beta-equilibrium gives an EOS for the core. In the inner crust, the NM EOS can be used 
to describe the dripped neutrons {x = 0) and the bulk matter in the nuclear clusters. A 
consistent model for the EOS of crust and core necessarily uses a unique NM EOS, and one 
should expect parameters characteristic of a given NM EOS to correlate with both crust and 
core properties. 

Nuclear matter with equal numbers of neutrons and protons {x = 0.5) is referred to as 
symmetric nuclear matter (SNM); nuclear matter with x = 0.0 is naturally referred to as 
pure neutron matter (PNM). Nuclei on Earth contain closely symmetric nuclear matter at 
densities close to nuclear saturation density />o ^ 2.7 x 10^^ g cm~^ = 0.16 fm~^ = no, 
where we use n to refer to baryon number density. Thus experiment has constrained the 
properties of E{r^ no, ~ 0.5) to within relatively tight ranges, but the properties of PNM 
remain uncertain from an experimental standpoint. In the past decade, much experimen- 
tal activity has been devoted to extending our knowledge of nuclear interactions to more 
neutron-rich systems and to higher and lower densities. Although we cannot produce pure 
neutron matter in the laboratory, we can produce matter with proton fractions as low as 
X ^ 0.3 in certain neutron rich isotopes and in the products of heavy ion collisions. This 
allows us to obtain information on how no, x) changes as x decreases. 

By expanding E{n^ x) about x = 0.5 using the isospin asymmetry variable 5 = 1 — 2x, 
we can define a useful quantity called the symmetry energy S(n), 



2 05^ 



(1) 



which encodes the change in the energy per particle of NM as one moves away from isospin 
symmetry. This allows extrapolation to the highly isospin asymmetric conditions in neutron 
stars. The simplest such extrapolation, referred to as the parabolic approximation (PA), 
truncates the expansion to second order, giving 



Ep^uin) = E{n, 6 = 1)^ Eo{n) + S{n) (2) 

i-np 
3no 



for the PNM EOS. Expanding the symmetry energy about x = where x = we 



obtain 

^(n) = J + Lx+^i^symx' + ..., (3) 

where J, L and i^sym are the symmetry energy, its slope and its curvature at saturation 
density. 

Since neutron star matter contains a low fraction of protons, many inner crust and global 
stellar properties are sensitive to the symmetry energy parameters J,L, etc. To give a sim- 
ple example, the pressure of PNM at saturation density is given by Ppnm(^o)=^o^/3. The 
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pressure in the inner crust and outer core is dominated by neutron pressure so a strong cor- 
relation exists between the pressure in neutron stars near saturation density and L. Neutron 
star EOSs which have higher pressures are often referred to as 'stiff; lower pressure EOSs 
are referred to as 'soft'. Thus, around 1 to 2no, 'stiff EOSs are associated with high values 
of L and 'soft' EOSs with low values of L. 

That large uncertainties exist in characteristic NM parameters such as J and L is one 
reason why many model predictions of potential neutron star observables span such wide 
ranges. On the other hand, observations of neutron stars offer the opportunity to obtain 
constraints on NM parameters and hence on the underlying models of the nucleon-nucleon 
interactions [11 J. 

The following is a non-exhaustive list of (potentially) observable neutron star phenom- 
ena whose precise properties depend on the properties of the inner crust and mantle and of 
the star as a whole; some of these will be described in detail elsewhere in this book. 

• Pulsar glitches. Young pulsars spin down under the action of magnetic torque, their 
rotational energy powering the radiation beam. Many are observed to undergo occasional, 
sudden, spin-ups called glitches [12J. Proposed mechanisms include crust-cracking as the 
star attempts to adjust its shape to become more spherical ifTSl and angular momentum 
transfer from one internal component to another such as some part of the crust superfluid 
neutrons to the rigid part of the crust |[T4l [HI, or a combination of both [|16i . In such 
models, the size, frequency and post-glitch relaxation of the spin period depend on, among 
other microscopic properties, the crust and core sizes, moments of inertia and composition. 

• Free precession. Certain pulsars exhibit long timescale periodic variation in their 
timing residuals, with periods of order years, suggestive of free precession of the star 
(TTKlSl. Free precession is can arise from mechanically or magnetically supported crustal 
deformation [20l |2l|, and the period depends also on details of crust-core coupling, 
notably through the properties of the crust and core superfluid l[22l[23l[24ll25l . 

• QPOs from SGR giant flares. Quasi-periodic oscillations (QPOs) in the tails 
of light curves of giant flares from soft gamma-ray repeaters (SGRs) have been observed 
t26||27l[28l[29|, and their frequencies lie in the range of possible torsional vibrations of the 
crust. The crust thickness, composition (through, e.g., the shear modulus) and the stellar 
size all affect the frequencies of such modes ll30l[3n[32l 13311341 . 

• Neutron star cooling. The crust thermalization timescale depends on the crust 
thickness as well as the thermal conductivity and specific heats arising from the heat trans- 
port mechanisms operating in the crust |[35l|36l. One intriguing possibility is the operation 
of the direct Urea process in the bubble phases of the mantle, where the delocalization of the 
protons may allow it [4]. The thickness of this layer plays an important role in determining 
how effective a cooling mechanism this might be. 

• Gravitational waves (GWs) from neutron stars. A rich array of stellar oscillation 
modes are possible, some of which might generate GWs detectable on Earth I^Ell, or 
lead to other observational signatures such as limiting neutron star spin-up |[39l|40|. Sta- 
bility of modes can depend sensitively on the physics at the crust core interface and crust 
thickness, ||4D|42|43l|44l|45l|43|43 GWs can also be generated by a quadrupole deforma- 
tion in the stellar shape supported, among other possibilities, by the elastic crust Il48ll49ll . 
Whether the crust is strong enough to support a large enough deformation to produce de- 
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tectable gravitational waves depends on the shear modulus throughout the crust, and thus 
its composition (especially in the mantle where the mechanical properties are particularly 
uncertain) and the crust thickness and stellar size ll33l . 

In this chapter we will review the dependences of the composition, thickness of the 
crust and the mantle, and certain global neutron star properties on the symmetry energy pa- 
rameters J and L. The interplay of such relationships in modeling neutron star observables 
and the potential for obtaining astrophysical constraints on J and L will be illustrated by 
constructing consistent crust and core models based on a model of uniform nuclear matter 
whose symmetry energy parameters can be smoothly varied. We shall use the compress- 
ible liquid drop model (CLDM) of crustal matter for its expediency; an outline of more 
sophisticated models will be given in the discussion at the end of the chapter. 



2. Symmetry energy and its correlations with observables 

In this section we review experimental and theoretical constraints on J and L and some of 
the correlations that emerge between J and L and neutron star observables. Throughout 
this chapter we will illustrate the correlations using a particular model for nuclear matter: 
the phenomenological modified Skyrme-like (MSL) model (501 [5T|. The energy as a func- 
tion of density n and isospin asymmetry 5 is written down in a form that closely resembles 
the uniform nuclear matter Skyrme EOS under the Hartree-Fock approach. The advantage 
of the MSL function is that its free parameters can be more easily related to the proper- 
ties of nuclear matter at saturation density and allows for a smooth variation of J and L 
independently while keeping the SNM EOS constant. The MSL EOS is written 

n / n \ ")^sym 

?loc /_\ _ /I ^.\-r^loc \ I ^.-r^loc / \ / \ 



where 



K7min) = (1 - y)Ell'M^^ + vK'rnino) i^-^J , (5) 



= — + n(Ceff + D^^5)', — - = — + n(Ceff - D^^5), (6) 



2m* 2m 2m* 2m 

77 = (3/5)(37r^)^/^, no is the saturation density and a, /3, a, Cgflp , D^^, 7sym, ^^syml^o) and 
y are free parameters. Ceff and i^eflf are fixed by setting the effective masses at saturation 
to be m* = 0.8m and m*o = 0.7m, where m is the average nucleon mass in free 
space. We set 7sym = 4/3; a, /3, a are set by the incompressibility Kq, density no and 
energy per particle of SNM at saturation, while E^^^^ and y control the absolute value 
of the symmetry energy and its slope respectively. We will keep constant no = 0.16 fm~^, 
E{n^, X = 0.5) = -16 MeV and Kq = 240 MeV. The symmetry energy in the MSL model 
is not restricted by the parabolic approximation (Eq. (2)), but includes contributions of 
orders higher than 5^ from only the kinetic and effective mass parts; the potential part of 
the symmetry energy, £^^^^(n), is quadratic in isospin asymmetry. 
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Figure 2. Left panel: Symmetry energy slope at saturation density L versus its magnitude 
Esymipo) = J extracted from isospin diffusion and neutron/proton (n/p) ratios of pre- 
equilibrium nucleon emissions within the Improved Molecular Dynamics (ImQMD-2009) 
model [|52l [53l [54l , isospin diffusion within the Isospin-Dependent Boltzmann-Uehling- 
Ulenbeck model (IBUU04- 2005) |[55l|56l, isoscaling |57J, energy shift of isobaric analogue 
states within the liquid drop model (IAS+LDM-2009) ll58]| , neutron(n)- skins of several 
heavy nuclei using the droplet model (DM) or the Skyrme-Hartree-Fock (SHF) approach 
El |59l [601, pygmy dipole resonances (PDR) in 68Ni,132Sn and 208Pb, f6Tl[S2l, and 
nucleon global optical potentials (GOP) |63 1. Taken from ref. |63|. Right panel: Correla- 
tions between J and L arising from the MSL EOS constrained by J =35 MeV, L=25 MeV, 
low density PNM calculations (together outlining our 'baseline' region) and ^(O.lfm"^), 
fits to nuclear masses using the Weiszacker |64 1 and Thomas-Fermi (TF) ll65l models, the 
Hugenholtz-van-Hove (HVH) theorem applied to global optical potentials f63\ from which 
the overall bounds on the left plot were extracted, PNM calculations taking into account 
uncertainties in 3-nucleon (3N) interactions using chiral perturbation theory ll66l and quan- 
tum Monte-Carlo (QMC) simulations |[67l|68l, the relativistic mean field (RMF) FSUGold 
model constrained by S'(0.12fm"^)=26MeV ||69l|7pl|, the RMF DD-PCl model constrained 
by S'(0.12fm-^)=27.6MeV EB, a best fit to a large number of SHF, Bruekner-Hartree- 
Fock (BHF) and RMF models Ell, and the model correlations from the ImQMD-2009 
analysis of isospin diffusion and n/p ratios from which the overall bounds on the left plot 
were extracted 1.53 J . Adapted from L73J . 

2.1. Experimental constraints on the symmetry energy 

The current status of our uncertainty in the symmetry energy based on a wide range of 
experimental results is summarized in the left panel of Fig. 2; the caption sunmiarizes the 
details. We point out again that the parameters J and L describe an idealized system, and 
the extracted values are somewhat model dependent. However, by using many independent 
observables we can check the consistency of our results and arrive at a more robust set 
of constraints. We will take as a conservative range from experiment 25<L<115 MeV, 
noting however that a convergence of experimental results to the lower half of this range 
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Figure 3. Pure neutron matter (PNM) energy per particle versus baryon number density 
from effective field theory (EFT) at low densities (SP) |79|, chiral EFT with 3-nucleon 
(3N) interations (HS) |[66l and quantum Monte-Carlo (QMC) with 3N interactions (OCR) 
16711681. The MSL EOS is plotted with symmetry energy slopes of L = 25, 70, 115 MeV 
and with J = 35 MeV CJ35'), ^(O.lfm-^) = 26 MeV CSq.i') (left panel) and constrained 
to fit the low-density PNM EOS of SP, OCR and HS (TNM') (right panel). Adapted from 
El. 

is apparent in Fig. 2. For J we take 25< J<35 MeV, which includes the ranges extracted 
from nuclear mass model fits 116411741 [7511761 . 

2.2. Theoretical constraints on the symmetry energy 

Symmetry energy constraints follow from theoretical calculations of the PNM EOS. At 
low densities (< 0.02fm~^) where only two-nucleon (2N) interactions need be considered, 
quantum Monte-Carlo (QMC), Green's function Monte-Carlo, effective field theory (EFT), 
chiral EFT and variational chain summation techniques 1 66} |67l |68l |771 |78l [79] have pro- 
duced robust constraints on the PNM EOS which may be taken as a constraint for our MSL 
model. In Fig. 3 we show the PNM constraint from Schwenk and Pethick (SP) [79 | as 
the red box at low densities; results since have converged on the lower bound of that box. 
Calculations at higher sub-saturation densities which require 3-nucleon (3N) interactions 
to be included start to diverge in their predictions; two recent calculations with estimated 
theoretical error bars for the 3N forces have been performed ||66l[68l (hereafter labelled HS, 
GCR respectively), and are indicated by the shaded bands in Fig. 3. 

2.2.1. Correlations between J and L 

In the left panel of Fig. 3, we plot the PNM EOS from our MSL model for J = 35 MeV 
and L = 25, 70, 115 MeV. A fixed J corresponds closely to a fixed value of the PNM EOS 
E{no, X = 0.0) at saturation density for a fixed SNM EOS (see Eq. (2)). Only the L = 70 
MeV EOS agrees with the predictions of theory at low densities. However, one can adjust 
J for a given value of L to obtain agreement with the low density PNM EOS; the results 
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for the MSL EOS are displayed in the right panel of Fig. 3. Doing so naturally introduces 
correlation between J and L; in the right panel of Fig. 2 we display the correlation obtained 
in this way for the MSL model. It is fit by J = 20.53 + 0.207L. For reference, the correla- 
tions obtained directly from the PNM calculations of HS and OCR, using the PA (Eq. (2)) 
with £^0 = —16 MeV to obtain J from £^pnm(^o), are depicted in Fig. 3; although offset 
slightly from the MSL results, their slopes are similar. A similar correlation is obtained 
from the Hugenholtz- Van-Hove (HVH) theorem which predicts a relation between J and L 
whose uncertainty can be related to global nucleon optical potentials ll63l 

One experimental probe of the symmetry energy is the measurement of neutron skins 
of nuclei. This probes the symmetry energy at densities around n = O.lfm"^; thus many 
models fix the symmetry energy at this density. In the right panel of Fig. 3 we show the 
MSL PNM EOSs constrained by 5(0. Ifm"^) = 26 MeV; varying L then produces a steeper 
correlation with J, also shown in the right panel of Fig. 2\ J — 29.0 + O.IL. It is worth 
noting that increasing the density at which one fixes the symmetry energy in a given model, 
increases the slope in the J-L plane. 

Similar correlations are obtained from two relativistic mean field models fTOl [Til and 
from a best fit to a wide selection of model predictions of J and L [72], also shown in the 
left panel of Fig. 2. Finally we also show correlations that emerge from nuclear mass fits 
ll64l|65]| and analysis of data from heavy ion collisions |53|. 

In what follows we shall use sequences of MSL EOSs generated by varying L with a 
variety of constraints on J: the sequence generated keeping J fixed will be labelled, e.g., 
'J35'; the sequence generated by fixing the low density PNM EOS will be labelled the 
'PNM' sequence; and the sequence generated by fixing ^(O.lfm"^) = 26 MeV will be 
labelled the 'Sq.i' sequence. The model correlations in the right panel of Fig. 2 overlap in 
the region 25<L<70 MeV, in line with the most recent experimental results. By combining 
the MSL 'PNM' constraint with the requirement that 25< J<35 MeV and L>25MeV we 
obtain a region in the J-L plane which we shall refer to as our 'baseline' region. 

2.3. Correlations with neutron star properties 

Some useful correlations of symmetry energy parameters with basic neutron star properties 
have been established, which we review here; more details can be found in the following 
references: [|IIl[70l|72l[80l[8ll 

• The pressure of neutron star matter in beta-equilibrium at no including the electron 
contribution can be approximated lITTllSTIl 



where the second term provides a correction of only 2-3% for L = 25 MeV, rising to 10- 
20% for L = 115 MeV, with J over the range 25 - 35 MeV. At densities slightly above or 
below this, extra terms are introduced, but the leading order will remain the one proportional 
to L alone. 

• The radius of a neutron star is found to correlate with the pressure at a fiducial density 




(7) 
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inside the star ifTTIl of the form 

i?^C(M,n)P(n)°-2^ (8) 

where the pressure in MeV fm~^ is taken at a density between no and 2no. C(M, n) is a 
constant for a given stellar mass M and fiducial density n. Given Eq. ^ to leading order 
and taking n = no, this suggests [821 

i??^0.5C(M,no)L°•^^ (9) 

using (no/3)^-^^ ^ 0.5. For a IAMq star, C(M, no) ^ 9.3. The compactness parameter 
^ = GM/Rc^ thus goes as 

6(M, nojc^ 

• Total moment of inertia of the star / scales with MB? dUSll, so 

I MC{M,no)L^-\ (11) 

• The crust-core transition density rice has been established to correlate inversely with 
the slope of the symmetry energy L72i[83l[84l[85l; if all other nuclear matter parameters are 
held constant, the correlation is generally found to be weakly parabolic: 

^cc = - br,L + CnL^, (12) 

where an - 0.1fm"^ 6n ^ 10-^fm"^MeV"\ Cn - 4 • 10-^fm"^MeV"^ depend on other 
model parameters for their precise values fEE\. For L ^10 MeV, the quadratic term is 
negligible, and the relation is closely linear; it becomes important for L ^100 MeV. 

• The inner crust-mantle transition density rip, that is, the density at which pasta 
appears, is roughly constant with respect to the symmetry energy parameters |[83]| . with the 
exact value depending on the inner crust model. Typically, rip ^ 0.04 — 0.06 fm~^. 

• The crust-core transition pressure Pec was initially though to correlate with L in a 
similar way to rice, as expected from Eq. Q; however, a more thorough survey of nuclear 
matter models reveal no significant correlation with Lor J [i70l|72l. This can be understood 
by the fact that although the pressure at a given density correlates positively with L, the 
transition density correlates negatively, and the convolution of the two correlations for a 
given nuclear matter model could be positive or negative depending on the other nuclear 
matter parameters. Over a wide range of NM models, therefore, there will be no obvious 
correlation. However, [72] did find a robust correlation relating Pec to symmetry energy 
parameters at densities characteristic of rice'- 

Pec = ap + 6p[L(0.1fm-3) - 0.343Ksym(0.1fm-3)], (13) 

where ap ^ 0.5MeVfm~^, 6p ^ O.Olfm"^ depend on other model parameters for their 
precise values ||72l[86l. Note also that Pec correlates linearly with the baryon chemical 
potential at the crust-core transition /Xcc 

• The relative crust thickness, mass and moment of inertia go as flmTOllSTIl 

AR R C 0.25,, .AM i?2 C2 . AI R' C\ 
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Higher order corrections to global crust and star properties depending on radius will be 
ordered in ascending powers of ^5; higher order corrections to global crust properties will 
be additionally ordered in ascending powers of Pcc/^ cc or equivalently jJjqq j ttIh where 
is the mass-energy density at the crust-core transition and is the neutron rest mass. 



3. The Compressible Liquid Drop Model 

In order to examine the symmetry energy dependence of inner crustal properties, we will 
use the simple, but physically insightful, compressible liquid drop model (CLDM) for the 
energy density of crustal matter. The CLDM formalism which we highlight below was 
originally formulated by BBP IS and updated to incorporate pasta shapes 1871 [88l. Its 
computational expediency makes it widely used as a model for the crust EOS (e.g. |[T0l[89l ) 
and suitable for a calculation of many crust compositions and EOSs over a range of nuclear 
matter parameters. Its simplicity is also its disadvantage, neglecting as it does important 
microphysics; these problems will be discussed later in the chapter. For further details of 
the model see ll2l 173118811. 

A unit cell of crust matter is approximated by an equally volumed cell with the geometry 
of the nuclear cluster under scrutiny; this is the Wigner-Seitz (WS) approximation. The 
dimensionality of the shapes is specified by a parameter d = 3, 2, 1 for spherical, cylindrical 
and planar geometries respectively. The total energy density of the matter an be written 
(neglecting rest masses) 

e^celK^c, X, n, rin) = V \nE{n, x)+5ex+£^th] +'?^feurf+£^(C+L)) + (l-'^)^n£^(^n, 0)+5e(^e), 

(14) 

where Tc is the radius of the WS cell, u = {r^jr^^^ is the volume fraction occupied by the 
nuclei or the bubbles of radii (or half- width in the case of slabs) tn, x and n are the proton 
fraction and baryon density of the charged nuclear component and the baryon density 
of the neutron fluid. is the number density of electrons. Charge neutrality demands 
rie = vnx where v is the volume fraction of the charged nuclear component, defined as 

^ _ { u nuclei 
[ 1 — bubbles, 

and the global baryon number density is related to the local baryon densities through rib = 
vn^ — v)n^. The contributions to the energy density of the cell break down as follows: 

• The electron kinetic energy density, that of an ultra-relativistic free Fermi gas, is 
= (3/4)/icfcene with fee = (Syr^ne)^/^. 

• The electrostatic contributions include the nuclear and lattice Coulomb energy den- 
sities, collectively written within the WS approximation as 

^(C+L) = 27T{exnrN)'^fd{u); fd{u) = 

The corrections to the nuclear Coulomb energy density from the finite surface thickness and 
the proton Coulomb exchange energy are given by 

4 3 

e^^{k, x) = --Tie^w^x^k^n; Sex{k, x) = - — 2^/^e^x^/^A:n, (17) 

9 4:7T 



d-2 



(16) 
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where k = (l.Syr^n)^/^ and it; is a distance representing the surface thickness, taken to be 
It; ^ 0.75 fm. Electron screening can be included but is neglected in what follows. 

• The energy per particle of the neutron fluid and of neutrons and protons in the bulk of 
the nuclear clusters, E{n^ x) and E{n^^ 0) are obtained using a model for uniform nuclear 
matter. We will use the MSL model outlined in the previous section. 

• The surface energy density 6surf can be written in terms of the surface and curvature 
tensions cTg, (Jq as 

^surf = rfcTs/^N + d{d - l)crc/rN, (18) 

Thermodynamic equilibrium constrains the surface and curvature tensions to be a function 
of only one free parameter, conveniently taken to be x |[90l|9Tl. One possible functional 
form is [HnHHlEl 

crs(x) = CTQ-^ — ; crc(x) = a^^a[l3 - x) + 7X . (19) 

XTP {1—X)P 

The parameters ctq, ctq^c, ^, o:, /3, 7 are typically adjusted to fit nuclear masses ESUH or 
microscopic calculations of the energy of the interface between two phases of semi-infinite 
nuclear matter (921 [93l using the same nuclear model responsible for the bulk nuclear matter 
terms in the CLDM. Finally, a description of a neutron skin can be consistently included in 
the CLDM 1 9111921 [9311941. 

The composition of the cell is obtained by minimizing the energy density of the unit 
cell with respect to the free parameters (e.g. x, n, rin and rc). This produces four equations 
to be solved which correspond physically to mechanical, chemical and beta equilibrium of 
the cell plus the nuclear virial relation ^(c+l) = 2esurf + ^curv which expresses the scaling 
between Coulomb, surface and curvature energy densities under equilibrium with respect 
to variation of the volume fraction of the charged nuclear component ||2|. 



4. Crust-core and inner crust-mantle transition 

We will now employ the CLDM formalism outlined in the previous section to calculate 
sequences of crust compositions and crust and core EOSs using the MSL nuclear matter 
model with L varied between 25 and 115 MeV. The crust EOS sequences will be labelled 
according the corresponding nuclear matter EOS sequences outlined in Section 2.2.1: we 
shall use 'J25', 'J30', 'J35',TNM' and Sq.i sequences. Results corresponding to the base- 
line region (Fig. 2) in J-L space will appear shaded in the plots that follow. 

The surface and curvature energy parameters in Eqs. [19] should be obtained either from 
nuclear mass fits or microscopic calculations for each different set of J, L values used. 
We adopt a different approach: we use the established correlation between the the surface 
symmetry energy (encoding the change in the surface energy in nuclei as they move away 
from isospin symmetry) and J (see, e.g., [80|) to fix h for a given J. Additionally we take 
(Jo = l.lMeV fm~^, p = 3 and (Jo,c = 0.6MeV fm~^ as typical values from microscopic 
calculations (see discussion in (731). The neutron skin is neglected. 

In Fig. 4 we show the variation of the crust-core and inner crust-mantle (= spherical 
nuclei-pasta) transition densities and pressures for the five sequences of MSL EOSs. For 
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Figure 4. Crust-core and spherical nuclei-pasta transition densities (left) and pressures 
(right) versus L for the 3 constant- J sequences 'J35', 'J30' and J25' and the sequences 
TNM' and 'Sq.i'. The effect of using the parabolic approximation (PA) to calculate crust- 
core transition properties for the TNM' sequence are also shown. The shaded region indi- 
cates the baseline results. Adapted from 11731 . 

the crust-core transition we also show results using the full parabolic approximation for the 
symmetry energy. We can summarize these results as follows: 

• The J25, J30, and J35 sequences. For constant J, the crust-core transition den- 
sity Ucc correlates negatively with L as expected from previous studies. Ucc varies by 
^ 0.06fm~^ in the range 25<L<115 MeV. For constant J, the crust-core transition pres- 
sure Pec correlates negatively with L, varying by ^ 0.6 MeV fm~^ over the same range, 
ncccorrelates positively with J, varying by ^ 0.04fm~^ over 25< J<35 MeV for constant 
L. A higher symmetry energy favors a larger proton fraction in the nuclear clusters, lower- 
ing their bulk binding energy more than their surface energy is raised, thereby making the 
clustered matter energetically favorable to higher densities. Pec also correlates positively 
with J for fixed L, varying by ^ 0.5 MeV fm~^ over 25< J<35 MeV. 

• The PNM and 5o.i sequences. When a linear, positive correlation of L with J, L( J), 
is imposed, the relationship of rice and Pec with L{J) convolves their relationships with L 
and J indepedently. Fig 4 demonstrates the evolution of the Ucc-L and Pqq-L trends as 
the slope of L{J) steepens from the 'J35' sequence through 'Sq.i' to TNM'; the negative 
slope of the correlations get less pronounced, and, in the case of Pec, becomes positive for 
the steeper L{J) slope of the TNM' sequence. 

• The parabolic approximation. The results for the PA diverge to higher transition 
densities as L increases compared to the full EOS. The MSL model is already parabolic in 
the potential part of the symmetry energy; our results thus demonstrate the importance of 
using the full kinetic part of the symmetry energy (as has been noted before 117211851 ). 

• Our baseline region gives a range of crust-core transition densities ncc=0.08 — 0.12 
fm"^ and pressures Pec ^ 0.35 - 0.85 MeV fm"^. 

• Relative to the crust-core transition density, the density of transition to the pasta 
phases rip shows little variation with L and J, tracing out a thin band np ^ 0.05 — 0.06 
fm~^; thus, the thickness of the pasta layers will correlate with L and J in a similar way to 
the crust-core transition density. The spherical-pasta transition pressure Pp increases with L 
for the TNM' sequence, whereas the variation with L at constant J is weaker. The baseline 
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Figure 5. Crust-core transition densities scaled by 30 MeV/J (left) versus L and transition 
pressures versus L{n = 0.1 fm~^) - 0.343 Ksy^{n = O.lfm"^) for the EOS sequences 
shown in Fig. 4. Adapted from 11731 . 

region gives Pp ^ 0.15 - 0.3 MeV fm~^. 

These results demonstrate that the crust core transition density and pressure are sensitive 
not only to the slope of the symmetry energy L at saturation but also to its magnitude J to 
a similar extent over the experimentally constrained ranges; thus, quite different results 
can be obtained from calculations using different NM models which adopt, implicitly or 
explicitly, different correlations between J and L. 

The transition densities for the sequence of constant J EoSs (Fig. 4) suggests a roughly 
inverse scaling of the densities with J. With this in mind, we plot in the left panel of Fig. 5 
the transition densities multiplied by 30MeV/J. The results form a relatively tight band 
correlated with L; included on the plot, as well as our MSL EOSs, are the results using 
a selection of Skyrme EOSs (see [73 1 for details). Similarly, using the suggested scaling 
given in Eq. (T3\ [72J, the right panel of Fig. 5 shows a tight correlation, reinforcing this 
result and indicating that more experimental data on the symmetry energy parameters and 
including the curvature Kgym would improve our estimate of the transition pressure, as well 
as hinting at the connection between the J-L correlations and the higher order symmetry 
parameters JCsym, etc. 



5. Crust composition 

We turn our attention to the composition of the inner crust/mantle versus density for a se- 
lection of EOSs from the TNM' and 'J35' sequences. Specifically, we show in Fig. 6 the 
volume fraction of the charged nuclear component v, the density fraction of dripped neu- 
trons Xn = (1 — i;)nn/nb, the local proton fraction of the charged nuclear component x, the 
WS cell size rc and two properties of crustal matter estimated from the compositional pa- 
rameters: the shear modulus /i ||95l|96l[97l at zero temperature and the melting temperature 
Tm (at which the crystalline lattice melts into a plasma) L98ll : 

/. = 0.1106 (f ) A-^/3</^l - X^r/HZef; = ^^g^. (20) 
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Figure 6. Top left: volume fraction v occupied by nuclei/nuclear clusters, top right: 
density fraction occupied by free neutrons, middle left: local proton fraction of clusters 
X, middle right: WS cell size rc, bottom left: shear modulus of crust scaled by pressure 
and bottom right: crustal melting temperature, versus baryon number density in the crust. 
Results are displayed for the L = 25, 70 and 115 MeV members of the 'J35' and TNM' 
sequences. The loci of the transition densities over the range L = 25 — 115 MeV are 
displayed for the TNM' sequence (larger black circles) and the 'J35' sequence (smaller 
blue circles). The results from the BBP [2] and DH |10| crust models are shown by the 
crosses and plusses respectively. The shaded region has the same meaning as in Fig. 4. The 
bottom two plots also show the results of the TNM' L=25 MeV EOS with a stiff surface 
energy. Adaped from from (731 . 

The volume fraction rises with density in the inner crust as the nuclei get larger and more 
closely spaced. The predictions of the various EOSs tend to converge close to the transition 
to the mantle, where the volume fraction approaches the ^ 0.125 estimate of the the Bohr- 
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Wheeler fission instability criterion (see, e.g. l[99|), shown here with the dashed line. Below 
this density, a higher L corresponds to lower volume fractions: the correspondingly higher 
dripped neutron pressure favors smaller nuclei. Above the mantle transition, predictions 
vary widely depending on L and the L — J correlation. The volume fraction reaches > 0.8 at 
the crust-core transition within the baseline region, a warning that the WS approximation is 
certainly not valid at the highest crustal densities. The crust-core transition fraction remains 
high for the TNM' sequence, but falls with increasing L for the 'J35' sequence down to 
0.1 (L=115MeV). 

The density fraction of dripped neutrons increases sharply with density at the top of the 
inner crust, before leveling off at intermediate densities to Xn ^ 0.7 — 0.9 with the higher 
values predictions of higher L. Close to the inner-crust-mantle transition, falls rapidly 
down to 0.02-0.2 at the crust-core transition, mirroring the behavior of v. 

The local proton fraction x generally decreases increasing density over the whole den- 
sity range. For J = 35 MeV, the crust-core transition fraction varies by about 0.03; from 
L = 25 MeV, x^c decreases from 0.06 down to 0.03 at L = 70 MeV, then increasing back 
up to 0.06 at L = 115 MeV; for PNM the variation is similar, starting at Xcc=0.025 for 
L = 25 MeV and increasing with L up to 0.06 at L = 115 MeV. In the lower density region 
the variation in x remains approximately constant at around 0.05. At a given density, higher 
J and higher L correlate with higher x. Higher J favors a higher proton fraction; higher L 
favors smaller, denser nuclei and nuclear clusters; as the synmietry energy increases with 
density, denser nuclei will also tend to favor higher x. 

The WS cell size rc decreases smoothly up to the transition to pasta, and then proceeds 
through a series of discontinuous jumps as matter transitions through the various nuclear 
shapes. In reality, these jumps may be smoothed out by the existence of intermediate shapes 
not considered in this work 1 100 J . Higher L, (higher neutron pressure) leads to smaller 
nuclei and smaller values of vq. 

For comparison, we also show the results of two of the most widely used CLDM crust 
EOSs, from Baym, Bethe and Pethick (BBP) [2J and Douchin and Haensel (DH) 1 10|. BBP 
over-predicts the strength of the surface energy of nuclei |101|; DH calculates the surface 
energy consistently with the bulk energy using the SLy4 Skyrme parameterization (L=45.5 
MeV, J=31.8 MeV). The DH EOS contains no pasta, a result of a slightly stiffer surface 
energy than in our baseline models (see discussion in 11731 ). 

The shear modulus is shown scaled to the crustal pressure and the melting temperature 
is given in terms of 0.01 MeV (1.16 x 10^ K). The baseline results for the shear modulus 
show a variation of a factor of about 2 throughout the inner crust; the melting temperature 
varies by a factor of 3 throughout the crust. The BBP and DH predictions are shown for 
comparison. The results from the L=25 MeV member of the 'PNM' sequence is also shown 
with a stiffer surface energy, which elevates the melting temperature by up to a factor of 2 
at high densities, but doesn't appreciably affect the baseline range for the shear modulus. 

6. Global crust properties 

An estimate of the mass and moment of inertia fractions of the mantle and the component 
of mantle containing the bubble phases (where the protons are delocalized) relative to the 
total crust amounts can be obtained using AMi/AMcrust ^ 1 — Pi/ Pec (see |[89l ) where 
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AMi is the mass and Pi the pressure at 
the upper boundary of the crustal compo- 
nent i = pasta, bubble phases. The ex- 
pressions for the moment of inertia frac- 
tions are identical. The mass fractions are 
shown in Fig. 7. For L < 70 MeV, the pasta 
phases account for between 50% and 70% 
the mass and moment of inertia of the crust, 
a result of those phases existing in the high- 
est density region of the crust. For L > 70 
MeV, the fractions remain above 40% for 
the TNM' sequence but fall off rapidly to 
just a few percent at L = 115 MeV for the 
'J35' sequence. 

The equivalent fractions for the bubble 
phases can also be significant: 10% - 25% 
for L < 70 MeV, remaining above 10 % for 
the TNM' sequence and dropping to zero 
at L 85 MeV for 'J35'. This region could therefore have an important influence on a 
range of phenomena, e.g. allowing the direct Urea process [4J. 

We extend our CLDM crust equations of state to the core using the same MSL EOS to 
calculate the pressure and energy density of beta-equilibrated nuclear matter in the core. At 
the highest densities, where a description of matter in terms of nucleonic degrees of freedom 
is expected to break down, we smoothly join the MSL EOS to two poly tropic EOSs of the 
form P = Ke(^+^/^) in a similar way to [ 102 J . This also ensures that our complete EOS is 
always sufficiently stiff at the highest densities to produce 2Mq neutron stars as demanded 
by observations II103II . The joins are made at energy densities of 300 MeV fm~^ and 600 
MeV fm~^ by adjusting the constant K to keep the pressure continuous at the join. The 
lower density polytrope has an index set at n = 0.5, while the second index takes a range 
n = 0.5 — 1.5 for values of L from 25-115 MeV respectively BTll . Note that although 
the additions of the polytropes substantially changes the maximum neutron star mass for 
small (soft) values of L, it does not substantially affect the radius and crust thickness of a 
neutron star of a given mass. Then, using the transition densities and pressures calculated in 
the CLDM model we can solve the general relativistic of hydrodynamic equilibrium (TOV) 
equations to obtain model static neutron stars and examine the global crust properties. 

Fig. 8 shows the percentage thickness of the whole crust and mantle and the crustal mo- 
ment of inertia compared to the equivalent global quantities. The trends of these quantities 
with L very much depend on the J — L correlation: for the TNM' sequence, the thicknesses 
and crustal moments of inertia rise monotonically with L, whereas for the 'J35' sequence, 
the relationships are distinguished by a rise with L up to a maximum at L ^ 45 — 70 
MeV depending on mass, and then a shallower decline with L beyond the maximum. The 
thicknesses change by a factor of ^ 2 over the range of L and by an order of magnitude 
from 1.0 - 2.0 M©; the mantle thickness is consistently about an order of magnitude lower 
than the total crust thickness. The relative crustal moment of inertia AI/I is a quantity 
of relevance in, e.g., the study of pulsar glitches. Under the assumption that glitches are 
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Figure 7. Fraction of crustal mass occupied 
by the mantle and the bubble phases therein 
for the TNM' sequence and three constant- J 
sequences of EOSs, with the baseline region 
shaded. Taken from f73\ . 
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Figure 8. Total crust and mantle thickness (left) and crustal moment of inertia (right) as 
a percentage of the total stellar radius and moment of inertia respectively. Results for the 
TNM' and 'J35' sequences are shown versus L for 3 stellar masses: 1.0, 1.4 and 2.0 Mq. 
Taken from |[86ll . 

self-regulating phenomena involving a constant time-averaged angular momentum transfer 
between some crustal component such as the dripped superfluid neutrons and the rest of the 
star, limits can be set on the minimum value of A/// for a particular pulsar given sufficient 
observational data ||l2l[T5 1. Current observational data on the Vela pulsar provides a limit of 
A/// > 1.6% im, which is shown as the horizontal dashed line on the right plot. Taking 
the above assumptions at face value, one can deduce from the TNM' sequence that stiffer 
EOSs are consistent with the Vela data for a wider range of masses; the softest EOS L=25 
MeV is consistent only for M <\25Mq. However, for the sequence 'J35', the L>85 MeV 
EOSs are inconsistent with the Vela constraint as are the L<50 MeV EOSs for M ^UMq. 
It should be emphasized that in determining the crustal thickness and moment of inertia, L 
alone is insufficient; they depend sensitively on how L correlates with J. 



7. Dependence of observable quantities on symmetry energy 

Finally we show two examples of simple estimates of potential observables that incorporate 
crust composition, thickness and global stellar properties discussed in the previous sections. 

The left panel of Fig. 9 shows the fundamental frequency 1=2, n=0 and first overtone 
1=2, n=l from the spectrum of torsional shear oscillations in the crust for a IAMq star 
(n is the number of radial nodes and / the angular constant associated with the spherical 
harmonics Y{^) L33J. Observed values of the frequency of QPOs at 16, 26, 28, 30, 84, 92, 
150, 155 and 625 Hz ||26l |23 EH |23 are indicated by the dashed lines. The calculation 
uses the value of the shear modulus at either the crust-core boundary, assuming the mantle 
is an elastic solid, or the inner crust-mantle boundary, assuming the mantle is a liquid (an 
upper limit on the difference in mechanical properties of the pasta from the rest of the 
crust). The TNM' sequence of EOSs is used. Superfluid effects are estimated through 
the relative mesoscopic neutron effective mass m*/m 110411 : a value of 1 corresponds to 
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Figure 9. Left panel: The frequency of the fundamental torsional oscillation mode (top) 
and the first overtone (bottom) in the the crust for a IAMq star as a function of L. The 
circles and triangles show the frequency assuming the shear modulus in the pasta phases to 
be that of a Coulomb lattice; the squares and diamonds show the frequency assuming the 
shear modulus in the pasta phase to be zero; the triangles and diamonds take into account 
the entrainment of the superfluid neutrons by the nuclear clusters. The dashed lines show 
candidate frequencies for the fundamental modes; 18, 26, 28, 30Hz in the fundamental 
frequency range and 84, 92, 150, 155Hz in the range of the first overtone. Right panel: 
Estimate of the maximum quadrupole ellipticity of a IAMq neutron star (top) and the 
corresponding gravitational wave strain (bottom) as a function of L. The strain is calculated 
assuming a neutron star frequency of 300Hz and a distance of O.lkpc from Earth. The filled 
circles indicate the value taking the shear modulus in the pasta phases to be that of a solid 
Coulomb lattice; the empty circles indicate the value taking the shear modulus in the pasta 
phase to be zero. The dashed line in the bottom plot indicates the sensitivity of the most 
recent LIGO science run [38]. Taken from |[33l . 

no superfluid entrainment effects. The frequency generically decreases with increasing 
L. Ignoring superfluid effects and the effects of the pasta phases, the frequency matches 
observed QPOs from SGRs only at the lowest values of L. If the pasta phases are purely 
liquid, the frequency falls by a factor of 3, making it difficult to match the 28Hz frequency 
observed, and being consistent with the 18Hz observed frequency only at the lowest value 
of L. One can see that, accepting the model and the interpretations of the results at face 
value, observed fundamental frequencies only match the predictions for L < 70 MeV. 

The right panel of Fig. 1 1 shows the maximum quadrupole ellipticity of the crust of a 
IAMq star, normalized to the canonical value of 10~''(crbreak/10~^), as a function of L 
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for the same TNM' sequence of EOSs. The equivalent gravitational wave strain is plotted 
in the bottom half for an 'ideal' reference neutron star with frequency 300Hz and a dis- 
tance of O.lkpc. The convolution of the trends of global star and crust properties and crust 
composition with L manifests itself in the non-monotonic variation of ellipticity with L, 
having a maximum at L p:^ 50-60 MeV. The sensitivity of the most recent LIGO science 
run is indicated by the dashed line; for this ideal neutron star, the GW strain from crustal 
mountains predicted in this simple model is detectable for L ^ 40-100 MeV offering the 
chance to distinguish between the various possibilities of crust properties. 

8. Conclusions and future directions 

We have constructed sets of neutron star EOSs that consistently encompass the inner crust 
and core, and include the crust composition and transition densities, using the compress- 
ible liquid drop model for the crust. These have been used to demonstrate the effect of the 
magnitude J and slope L of the symmetry energy at nuclear saturation density on micro- 
scopic and global crustal properties, and potential neutron star observables. The crust-core 
transition density and pressure, crustal composition, stellar and crustal mass, thickness and 
moment of inertia, torsional crust oscillation frequencies and maximum crust deformation 
all depend sensitively on both J and L within their experimentally constrained ranges. One 
of the dominant neutron star model dependences is therefore the correlation between J and 
L which constrains their possible values in J — L space. Experimental and theoretical in- 
formation about L and J and their correlation will continue to improve; in order to add 
neutron star observations to this investigation, consistent explorations of neutron star prop- 
erties over the constrained ranges is of great importance. The sets of EOSs used in this 
paper are available to interested parties II105L 

The simplicity of the CLDM allows useful exploration of the dependence of compo- 
sition and transition density on J and L, but it possesses several drawbacks. Firstly shell 
effects within the nuclei, or arising from scattering of dripped neutrons off of nuclear clus- 
ters, are ignored. Such shell effects can dominate the determination of nuclear geometry, 
the equilibirium size of the nuclear clusters (which will proceed in discrete jumps corre- 
sponding to changes in the nuclear 'magic' numbers with density) and the ordering of the 
pasta phases, as well as transport properties such as contributions to heat transport from nu- 
clear components II 1061 and entrainment of dripped neutrons by clusters II104L Secondly, 
the WS approximation is expected to break down at when the nuclear separation becomes 
comparable to the cell size II107L which occurs in the mantle. Thirdly, effects that act over 
ranges greater than the unit cell are not consistently accommodated in the CLDM; longer 
range electron screening, larger scale self-organization of pasta phases and long wavelength 
transport effects are all unaccounted for. 

Some of these effects can be taken into account by more sophisticated crust mod- 
els. The Thomas-Fermi and Extended Thomas-Fermi methods (e.g 111081 11091 IllOII ) are 
semi-classical models employing the local density approximation which allows the nuclear 
surface energy to be calculated self-consistently with the bulk nuclear energy. Shell cor- 
rections can be self-consistently added using the Strutinsky integral method llll l. The 
ID-Hartree-Fock method is fully microscopic and accounts for surface and shell energies 
self-consistently (e.g. 111061 11121 11131 11141 ). but it describes only spherically symmet- 
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Figure 10. Integrated neutron density in a 
cubic cell at densities of n=0.04, 0.06, 0.08 
and 0. 1 fm~^ (top left to bottom right) cal- 
culated using the 3DHF method [73] with 
the SLy4 Skyrme parameterization. 



Figure 11. Transition densities obtained 
using the 3DHF method compared to those 
obtained using the dynamical method of 
t85| for the SLy4, SII, SkM* and SkMp 
Skyrme parameterizations. 



ric configurations and is thus constrained by the spherical Wigner-Seitz approximation. 
This restriction is lifted in the more computationally demanding 3D-Hartree-Fock (3DHF) 
method Illl5[|116i , allowing a self-consistent probe of the shape-phase- space of pasta lay- 
ers with shell effects included. As example, some results of a 3DHF model 111171 are shown 
in Figs. 10 and 11. The former shows the local neutron density is plotted over the unit 
cell at four different densities encompassing the crust-core transition, showing the evolu- 
tion of nuclear shape from spherical through bubble to uniform matter. The latter shows 
the resulting transition densities for four different parameterizations of the Skyrme inter- 
action, compared to the dynamical method of |85|. The Hartree-Fock method is naturally 
extended to include pairing effects self-consistently (Hartree-Fock-Bugoliubov). Longer 
range effects can be simulated via the semi-classical quantum molecular dynamics method 
Illl8[|119lll201ll21ll . Between these models, a complete physical description of the crust 
can be built up. Many of the methods mentioned are much more time-consuming than the 
CLDM, making a wide-ranging survey of nuclear matter parameters unwieldy. However, if 
we know the quantities for which the CLDM provides a reasonable estimate, and to what 
densities it remains reasonable, one can use the CLDM as a useful guide for more realis- 
tic calculations. Ultimately, the goal should be to have all relevant microscopic inputs to 
neutron star models calculated consistently with the nuclear matter EOS; much work still 
needs to be done in this direction. 
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